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ABSTRACT 


A  theory  of  electrolytes  founded  on  the  principles 
jf  statistical  mechanics  is  presented.     The  primary  physi- 
cal assumptions  are:     1)  point  ions;  and  2)  small  fluctuations. 
The  theory  is  a  modification  of  Kramers'    theory  of  electro- 
lytes.    The  modification  consists  of  a  more  physically  accu- 
rate development  of  assumption  2).     The  main  result  of  the 
present  theory  is  that  a  partition  function  for  the  electrolyte 
is  obtained  for  all  concentrations     whereas  Kramers'    theory 
breaks  down  beyond  a  limiting  concentration.      The  deviation 
from  the  Debye -Hueckel  limiting  law  is  in  the  direction  of;  the 
experimental  observations. 


INTRODUCTION 


In  the  early  days  of  the  development  of  the  theory  of  strong  electrolytjek, 
Kramers      showed  that  the  Debye -Hueckel  low  c  jncentration  limiting  laws  could 


1.     H.A.    Kramers,    Proc.    Roy.   Acad.   Amsterdam  XXX,    145  (1927). 


be  derived  from  the  statistical  mechanical  partition  function  as  well  as  from  the 
Po  iss  jn-Boltzmann  differential  equation.      However,    with  increasing  concentration 
the  Kramers  and  experimental  results  diverged  in  opposite  directions  from  the  Debye- 
Hueckel  limiting  laws.     This  can  be  seen  in  Fig.    1  where  the  osmotic  coefficient 
is  plotted^  as  a  function  of  the  square  root  of  the  concentration  for  a  (l-l)  salt,    a 
salt  of  NaCl  type. 


2.  The  osmotic    coefficient  was  defined  by  N.    Bjerrum  as  follows.     Let  P  =  NkT  be 
the  osmotic  pressure  of  a  solution  of  N  independent  (non-interacting)  ions.     Also,    let 
P  =  gNkT  be  the  osmotic  pressure  of  the  corresponding  real  solution  in  which  the  ions 
interact  with  each  other.     Here  g  is  the  correction  factor  to  P  for  a  real  solution. 
The  ratio  (P   -  P)/P  =   1   -  g  is  the  osmotic  coefficient. 

3.  Detailed  accounts  of  the  observed  behavior  of  electrolytes  and  of  the  consequen- 
ces of  the  Debye -Hueckel  theory  are  given  in  the  two  treatises:     a)  H.    Falkenhagen, 
Electrolytes,    Oxford  (1934);     b)  Harned  and  Owens,    The  Physical  Chemistry  of  Elec- 
trolytic Solutions       Reinhold  (1950).     Very  few  references  to  the  classical  papers 
w.ll  be  given  here.     They  can  be  found  in  a)  and  b). 


An  especially  interesting  feature  of  Kramers      theory  is  the  existence  of 
a  critical  concentration  above  which  the  partition  function  becomes  a  complex 
number.     To  understand  Kramers'    interpretation  of  this  result  let  us  suppose 


that  we  have  obtained  an  exact  formula  for  the  osmotic  coefficient  of  a  system  of 
ions  of  ionic  radius  a.     A  curve  such  as  one  of  those  plotted  inE  g.   I    could  be 
drawn  for  each  value  of  a.     In  particular  there  would  be  a  limit  curve  as  a — >■  0. 
If  the  result  of  Kramers  were  exact,    the  critical  concentration  would  be  that 
concentration  above  which  the  limit  curve  would  not  exist.     Physically  this  might 
mean  that  if  the  concentration  of  a  dilute  solution  of  point  ions  were  gradually  in- 
creased,   a  sudden  association  of  oppositely  charged  ions  would  occur  at  the 
critical  concentration  with  a  correspondingly  sudden  drop  in  the  osmotic  coeffi- 
cient. 

The  Kramers  critical  concentration  is  quite  low,    and  corresponds  to  0.03 
mol  per  liter  for  KC1  in  water  at  15    C  and  to  0.  0005  mol  per  liter  for  MgSCK, 
Even  though  real  ions  have  non-vanishing  ionic  radii,    one  would  expect  that  if 
the  association  divergence  in  the  osmotic  coefficient  for  point  ions  were  correct, 
then  some  residue  of  it  should  still  appear  in  the  experimental  curves  in  the 
neighborhood  of  the  critical  concentration.     Actually  the  experimental  deviations 
from  the  limiting  laws  for  some  "small  ion"  salts  (KNO,)  do  not  become  sig- 
nificant until  the  concentration  is  considerably  higher  than  the  critical  one. 

Since  Kramers'    paper  contains  several  mathematical  approximations, 
to  which  he  called  attention,    we  have  felt  that  it  would  be  of  interest  to  see 
whether  or  not  his  qualitative  results  depend  on  the  mathematical  approxima- 
tions.    By  replacing  what  we  consider  to  be  Kramers'    most  severe  approxima- 
tion by  one  that  is  physically  more  reasonable,    we  shall  show  that  the  diverg- 
ence at  the  critical  concentration  disappears  and  that  the  new  deviations  from 
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Fig.    1.      The  osmotic  coefficient,    1    -  Q   ,    as  a  function  of  the  square  root  of  the 
concentration,    C     for  a  (l-l)  electrolyte,     a)  The  Kramers'    Tneory:     Cj^  =   .029 
moles/liter,     b)  The  Modified  Kramers'    Theory:     CMT^  =   .78  moles/liter,    c) 
The  Debye-Hueckel  Limiting  Law.    d)  Scatchard  and  Prentiss     J.   Am,    Chem. 
Soc,    44,    2690(1932).    e)  Scatchard  and  Prentiss,    J.   Am.    Chem.    Soc,    55, 
4355  (1933). 
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Fig.    2.      The  osmotic  coefficient,    1    -   9  ,    as  a  function  of  the  square  root  of  the 
concentration,    C,  for  a  (2-2)  electrolyte,     a)  The  Kramers'    Theory:     Cj^  =   .00048 
moles/liter,     b)  The  Mo  dified  Kramers'    Theory:     CwK  =    ,013  moles/liter,    c) 
The  Debye  Hueckel  Limiting  Law.    d)  Hovorka  and  Rodebush,    J,    Am.    Chem.    Soc, 
49,    647  (1927). 


the  Debye-Hueckel  limiting  laws  (for  point  ions)  are  in  the  same  direction  as  the 
experimental  deviations.     Inasmuch  as  we  do  not  include  short  range  repulsive 
force  parameters  in  the  theory,    we  obtain  a  single  function  for  thermodynamic 
quantities  of  salts  of  a  given  ionic  type  rather  than  the  family  of  functions  that 
has  been  observed  experimentally. 

KRAMERS'    THEORY  OF  ELECTROLYTES 

We  shall  present  the  main  features  of  Kramers'    theory.     This  is  done 
because  Kramers1    really  elegant  analysis  is  not  widely  known  and  because  our 
modification  will  then  be  readily  understandable. 

Let  N  ions  with  charges    €,     £^  , ,€w   be  located  at  positions 

r   ,  .  .  .  ,    r-*,  in  a  container  of  volume  V.     The  thermodynamic  properties  of  this 
system  can  be  computed  from  the  partition  function 

v 

Here  df  .  is  the  volume  element  of  the  jth  ion  and  U("r .  ,  .  .  .  ,  rN)  is  the  total 
potential  energy  of  interaction  of  all  the  ions.     The  symbol  A\i)     denotes  the 
deviation  of  the  free  energy  of  the  system  of  ions  from  that  of  a  system  of  N 
independent  point  particles. 

It  was  discovered  by  early  workers  that  if  U(r    ,  .  .  .  ,  rN)  represents 
only  Coulomb  interactions, 


(D  =  dielectric  constant  of  the  solvent  and  r       is  the  distance  between  the  jth  and 

Jk 

kth  ions),    then  the  thermodynamic  quantities  computed  on  the  basis  of  (l)  di- 
verge.    This  short  range  divergence  does  not  really  occur  because  of  the  exist- 
ence of  short  range  repulsive  forces.      In  the  work  of  Kramers  the  divergence  is 

avoided  by  essentially  excluding  ionic  configurations  involving  very  saort  distances. 

4 
Kramers'    analysis  is  based  on  a  theorem  of  O.    Klein     which  states  that 


4.     O,    Klein,    Medd.    K.    Vet.    Nobelinstitute,    5,    No.    6,    (1919). 


if  U(r ,  ,  .  .  .  ,  rj  is  a  homogeneous  function  of  degree  |j.,    that  is, 
U(\  r .....  ,    \rjSr)  =  X.^U(r.     .  ,.,rj,    then,    for  N  sufficiently  large, 

e^O^/feT]    =  [F(r)-'y3AT)]N 

where  n  is  the  number  of  particles  per  unit  volume  and  F  is  some  function  of 
(n  /kT),     Klein'  s  Theorem  can  not  be  applied  if  the  ionic  radius  or  any  of 

the  usual  expressions  for  the  repulsive  forces  is  introduced  because  the  homo- 
geneous character  of  the  potential  energy  function  is  lost. 

The  Coulomb  interaction  as  a  function  of  distance  is  an  homogeneous 
function  of  degree   -1.     In  order  to  avoid  the  short  range  divergence  and  pre- 
serve the  homogeneity  of  (2),    Kramers  replaces  (Z)  by 


(3) 


where 


?ik  =  ft"'!]  -«K-AW]    > 


As  rjj^— >0,    g.,  — ^0,     However,    if  X.    is  chosen  to  be  a  very  large  number,    then 

§ik  =  rik        except  when  r^  is  extremely  small.     By  locating  the  origin  of  the 

1/3 
coordinate  system  outside  the  volume  V,    R  .,    is  of  order  V    '     .     Since  the  mean 

distance  between  ions  is  of  order  n~    '     ,    n  =  N/V,    h.,    is  of  order  N  .     Hence, 

by  choosing  \   to  be  of  order  N    '      the  potential  energy  function  (3)  does  not  dif- 
fer significantly  from  the  Coulomb  function  even  for  r  •,    considerably  smaller 
than  the  mean  distance  between  ions.      When  there  is  considerable  association  of 
oppositely  charged  ions,    i.e.   at  high  concentrations,    (3)  cannot  be  expected  to 
give  quantitative  results.     Since  the     constant  X.    does  not  appear  in  the  final  re- 
sults,   there  is  no  way  of  varying  the  thermodynamic  properties  with  the  char- 
acteristics of  the  repulsive  forces. 

In  his  detailed  calculations,    Kramers    does  not  use  (l)  but  a  partition 
function  analogous  to  the  Gibbs  grand  partition  function.     His  generalized  en- 
semble is  such  that  for  a  given  configuration  of  the  centers  of  the  N  ions,    the 
charge  on  a  given  ion  may  assume  any  value  which  the  charge  of  an  arbitrary 
ion  can  have.      The  new  partition  function  is  defined  by 

V^yf,[-^AT]=HJVA  N"N 2t*ftU/kT]  ,      <*> 

v 


where  the  summation  under  the  integral  extends  over  all  possible  arrangements 
of  the  charges  £    >£»,  "  "  "  )  ^N     on  the  i°ns  located  in  the  volume  elements 

cJr       dr  .   _     d?       5 


5.     An  investigation  of  thermodynamic  fluctuations  in  this  generalized  ensemble 
has  been  made  by  S.   Levine,    Proc.    Roy.    Soc.   A152,    529  (1935).  • 


The  total  potential  energy  is  a  quadratic  form  in  the  charges      j  £•  C    • 
Kramers  diagonalizes  this  quadratic  form  by  introducing  a  set  of  variables  }  y.  ( 
which  are  linear  combinations  of  the  variables     }€ :  (  s         ^l  =    2_    "fc  !  ^j    ^ 
and  choosing  the  matrix        f  )fw  •  J  so  that  U  =  |    2_-    b-y-    .      The  b'  s  are  the 

characteristic  values  of  the  matrix  (D''Ql»)  . 

In  order  to  take  advantage  of  the  diagonalizaf  ion,    Kramers  then  assumed 
that  the  variables  I  y-  y   could  be  regarded  as  continuous  variables  and  that  the 
distribution  of  the  set  I  y;  (    was  Gaussian  in  each  variable.     The  probability 
that  yi    lay  in  the  range  dy,    was  given  by 

(iTrlT'/2e4[-y^/2^]dyk    • 

Here    f)      is  a  number  proportional  to  the  ionic  strength  and  is  the  standard  de- 
viation from  neutrality  of  a  charge  chosen  at  random  in  a  solution. 

It  was  noted  by  Kramers  that  the  characteristic  values   [  b^  ,    which  de- 
pend on  the  ionic  configuration,    are  not  all  positive  for  every  configuration. 
This  implies  that  for  T  sufficiently  small  integration  over  y,  ,    say,    would  lead 
to  an  infinite  value.     It  is  the  non-positive  definite  quadratic  form  coupled  with 


the  assumption  of  statistically  independent  y^'  s  with  Gaussian  distribution  which 
is  responsible  for  the  breakdown  of  Kramers'  theory  above  a  critical  concentra- 
tion.    Increasing  concentration  is  equivalent  to  decreasing  temperature  because, 

from  Klein*  s  T  heorem,    these  two  variables  are  equivalent  to  the  single  variable 

1/3  At 

n        /T . 

This  concludes  our   resume  of  Kramers'    work.      Further  details  involve 
the  evaluation  of  the  partition  function  (4),    and  these  details  will  appear  below  in 
the  evaluation  of  our  modified  partition  function. 

THE  MODIFIED  KRAMERS'    THEORY 

An  important  parameter  of  a  system  of  ions  is  the  reciprocal  length    X, 
which  is  defined  by 

K2  =(41r/DfeTj{(N;/V)e2if  (5) 

where  s  is  the  number  of  ionic  species,    N^  is  the  number  of  ions  of  the  ith  species 
and  ez-  is  the  charge  on  an  ion  of  the  ith  species  (e  being  the  electron  charge). 
Let 

rL=DI?TVx,2/4T7M    •  <6> 

This  is  a  constant  proportional  to  the  ionic  strength,    and  we  may  write,    using 
equations  (5)  and  (6), 

N  2 

X  £a    =N^    .  (7) 


Our  modification  of  Kramers'    theory  consists  of  computing  thermody- 
namic averages  over  the  ensemble  restricted  by  Equation  (7)  for  fixed    n  , 
A  charge  configuration  in  our  ensemble  will  consist  of  a  set  of  N  ions  with 
the  properties: 


i/a  j/2 

0       -(Nr|)     ^   £^    (Mr|)  for  all  j; 


z 


b)       iGj-Nrj       ^ 


the  value  of  the  sum  it\  the  real  physical  system  of  interest.     Our  complete  en- 
semble will  be  the  set  of  all  systems  of  N  ions  with  the  proper  ionic  strength 
(i.e.    with  the  proper  value  of  Nfl  )  which  are  confined  to  a  container  of  fixed 
shape  and  volume  V  at  constant  temperature  T. 

The  average  value  of  a  function  of  position  and  charge,    F(  \  r  .  I  >  jjS;  j  J  i 
will  then  be 


<f>»Z"'M  l'W^A'^F^U^I>tu/^  » 


(8) 


where  Z  is  the  partition  function  given  by 


7  =vn«|=[-^t]=aS'"5  S-^v<Mv<^  «t>[-uM] .  (,) 

The  constant  A  is  chosen  so  that  as  T-*-OQ   .    AU*  -*"   O       .     A  is  the  recip- 

i_ 
rocal  of  the  surface  area  of  an  N-dimensional  sphere  of  radius  (NO)2. 


tf.zv^&^/rw)  . 


(10) 


Before  attempting  to  compute  the  above  partition  function  let  us  compare 
the  modified  ensemble  with  that  of  Kramers,      Kramers'    basic  partition  function. 
Eq.    (4),    may  be  considered  to  be  correct  (to  within  the  limitations  inherent  in 
the  manner  in  which  the  Coulomb  divergence  is  avoided  as  r  .,  ->- 0)  at  least  for 
symmetrical  electrolytes.        However,    after  making  his  mathematical  approxi- 


*See  Ref.    5,     Apparently  in  unsymmetrical  salts  some  corrections  might  have 
to  be  made  in  the  free  energies. 


mations,    Kramers'    ensemble  becomes  one  in  which  the   I  £  :  \    are  not  limited 
to  discrete  values  ezj  ,  .  .  ,  ,  ezs  but  are  allowed  to  range  between  -  qO    and  +  00 
with  the  probability  for     6:    to  lie  between      £•        and    £.     -j-  uGj       given  by 
(2TT0")      eX.bj-€j/2n  J  d€j     .  Averages  are  taken  over  all     £^S 

and  all  r,  '  s.     Since  the     1  £.  \      are  independent  variables,    the  ionic  strength 
for  the  systems  in  this  ensemble  ranges  from  0  to   QQ     .     The  modified  ensem- 
ble  is  that  part  of  Kramers"    ensemble  which  contains  systems  of  the  correct 

ionic  strength.      The   €:  $    are  not  independent  but  connected  through  the  rela- 

N       2  I 

tion         2_    £\     =    •'  0  .     It  is  this  property  which  makes  the  modified  en- 

semble  more  representative  of  a  real  electrolyte  than  Kramers'    ensemble. 

To  evaluate  the  partition  function,    Eq.    (9),    we  (as  did  Kramers)  take 
advantage  of  the  fact  that  U  is  a  quadratic  form  in  the  charges.      Multiplying 
the  integrand  Of  Eq     (9)    by  the  normalized  delta  function 


10 


H 


^.^.ifflal'Jh^fUl)]  .       (n» 


-00 

we  have 

00 


(12) 


Z^k\]  j--Jck1..(l€Ndr...drNSC€b...)Oft^[-u/fcT3  • 

-oo      v 

Rearranging  the  order  of  integration,    we  write 

Z  =  A  if'  (N  /  j dx  ex|{-  ix  N {| j  jj?,  J^  jk,-  <^  a$xl  f  -  «*I^J  e^_ 

where  the  line  ixn  =  constant,    xq  real  and  positive,    is  above  the  singularities 
of  the  integrand.     The  rearrangement  is  possible  because  U  as  defined  in  (3) 
is  finite  for  all  positional  configurations  of  the  ions.      Therefore  an  x-domain 
will  exist  for  which  the  quadratic  form  in  the    j  £.    (    is  negative  definite.      Thus, 
the  integration  over  the  )  £..  C   will  be  finite.     We  then  have  that 

lot)  L    i=l  i<k 

where   D^  is  the  determinant  of  the  quadratic  form  in  the  exponent. 
Therefore 


11 


lXo+00 


2= A(^/2^N"5d,<e^[-ixNri]j"5^V"dTM  •  w 

lX0-  00  V 

We  shall  briefly  digress  in  order  to  point  out  that  the  condition  of  elec- 
trical neutrality  holds  in  the  mean  as  it  does  in  Kramers'    theory.     The  total 

charge  of  a  system  is  Q  =/     C:  so  that   \0/    —  Ns^C.:/,        Now  a  linear, 

orthogonal  transformation  from  the  variables    )£•   (    to  variables  1  y-  I   exists 
such  that  the  quadratic  form  in  the  )  c:  (  becomes       £_  b.y-    ,     This  means 

J 

K 

that  \G>/  =    2_    *jk  \^k/    v      the  Y>L  S      being  the  transformation  coeffi- 
Fr  ™ ' 

cients.     Since  the  integrand  in  Eq.    (13)  becomes  an  even  function  of  y.  for  all 

j,    then  ^yj^  =  0  for  all  j  and  <Q>  =   0. 

Now  we  may  write  J  Y  +  oO 

ex|>[-AV>M>  A(Nrlf  TTiN"'  j  d*  «|.[-uNn]  GN  ,    (17) 

ixc,-"50 


with 


-•/2V  (18) 


V 

Following  Kramers,    if    £\*\Jj     is  to  have  thermodynamic  sense,    G^  would 
be  the  most  probable  value  of  Djsj        as  N  ->"  o0    >    where  x  would  be  considered  as 


12 


a  parameter  of  external  constraint  on  the  system.     Assuming  that  fluctuations 

are  small,    which  is  valid  for  not  too  high  concentrations,    Eq,    (18)  can  be 

written  as 

1 

GN=<Dn)  <19> 

V 

This  means  that  as  N  — >■  00     ,    G^  is  the  reciprocal  square  root  of  the  most 
probable  value  of  Dj^. 

If  we  followed  Kramers  in  applying  Klein1  s  Theorem  to  Eq.    (18),    we 
would  have  GN  =       I  CD  (x,    n         /T)  1       for  N  sufficiently  large,    where    <p     is 
some  function  of  the  variables  x  and  n    '     /T.     This  approximation  is  not  suf- 
ficiently exact  for  our  purposes  because  the  essential  singularity   of  Gj^  is  not 
properly  characterized  by  the  function  <T>    .     The  essential  singularity  of  G-^  can 
be  put  in  evidence  by  noting  that  the  transformation  from  the  variables    )  £.   ( 

{  y]   is  such  that        Z     ?j^j^  *  £    ^J  ^J      ' 
u  J^K  J  =  l 


to  the  variables 
We  can  then  write 


N       j- i 


It  is  plain  that  we  may  write 


N  - 

d  =  TT[-tx  -r(2DNTrlbJj 

N  •    ,     I— 


(20) 


<^X=   Tl[-ix  +  (2DhTrAj]  ,  (21) 

where  (2DkT)_1X.  •  is  one  of  the  N  roots  of  the  polynomial    \E>N  ^>    y.     Since 
the  matrix     (Q  -u>)  is  real  and  symmetric,    the  roots  are  real  and  \\  .  |    ^     X.  rn(N) 
where  \      (N)  is  the  magnitude  of  the  maximum  root  of  \Dj^  ^>   v-     The  maximum 
root  may  be  independent  of  N.     This  also  means  that  in  Eq.    (17)  we  must  take 
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x     N   (2DkT)_1\      (N), 
o  '  m 

Let  us  put       oC~      =   -ixZDkT,    and  write 

<.VN\=  <,-i*?ex\>\H%u]     >  (22) 

where  N 

The  function     Q    (oC)       is  analytic  in  the  region      |  C<  |   <C     1A      (N). 

Thermodynamic  arguments  require  the  existence  of     Q     (bO  as 

N*->-   Oo      .    that  is,     Q      (jpM     must  exist.     Since  X.      (N)  exists  for  every  N     and 

since  the  radius  of  convergence  is   1/X.      (N)  for  every  N,    the  existence  of  the 

m 

limit  function     Q.fo(j  implies  a  non-zero  radius  of  convergence  and  the  conse- 
quence that  X.      (N)  — >*  X.      (OQ)    =     X      ,    where  a  finite  X.        must  exist.     Therefore, 
m  m  m  m 

X.        may  be  determined  from  the  radius  of  convergence  of    Q^lpC)  if  this  func- 
tion can  be  found. 

The  existence  of  a  finite  X.        has  the  following  consequence.     Since 

m  51 

IX.  .  1^  XiN^  f°r  a^  J'    the  N  roots  are  confined  to  a  finite  real  interval.     In 
general,    the  N  roots  will  be  all  different.     Therefore,    as  N— ^-  OO      ■    the  roots 
will  become  dense  in  the  limit.     It  then  follows  that  the  behavior  of      9<y\  C°0 
as       (X_  -^-    1/X      ,    that  is,    the  character  of  the  singularity  of      GL-.     on  the  ra- 
dius of  convergence,    depends  on  the  limiting  density  of  the  X.;.     It  is  then  clear 
that  the  character  of  the  singularity  of     Q  p~    has  little  to  do  with  the  single 
term  In    1  1  +  0(.  X.      (N)  I     in  the  sum      Q..    which  exactly  describes  the  singula- 
rity at    0(_  =  Up    l/X.rn(N)  for  every  N       This  means  that  an  accurate  description 
of  N  Q     .        as  N   — >■  OO       requires  us  to  write 
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CN  (23) 

where  f\j 

?«>-ir.iir£,n(,+eaj)    > 

and  we  have  supposed  that  the  singularity  occurs  at     <A.  = -fl  /X.       .     Because  the 


m 


roots  become  dense  in  the  limit,      Q  (<?C)  —  9^  v^/»    Hence,    we  are  required  to 


write 

<"DN>  =(|-*Aj(-U)Nex|=[Ng]  (24) 

as  N -^    CX)       •     It  should  be  noted  that  the  previously  mentioned     Cp  =  exp  Q  . 
This  more  accurate  limiting  expression  for  <\DKT/  v  does  not  influence  the 
result  obtained  by  Kramers. 

The  application  of  Klein'  s  Theorem  implies  that    Q     will  be  a  function 

of  the  two  variables  x  and  n         /T,    or    (?^  and  n         ,    and  that  X.        will  be  pro- 

♦  ■        i  +       1/3 
portional  to  n  ' 

Eq,    (24)  is  equivalent  to  stating  that  for  N  large,      \Pn/  v  =  O(M^) 

where  M  is  a  function  independent  of  N.      We  followed  Kramers'    thermodynamic 

reasoning  to  obtain  this  result.     It  is  implicit  in  Klein5  s  Theorem  that  the  Nth 

root  of  the  partition  function  integral  be  bounded  as  N— >»  oQ  It  is  clear, 

however,    that  the  asymptotic  behavior  of  D^.,    and  that  of    \D     /  v,    depends  on 

the  properties  of  the  elements  of  the  determinant.     For  example,    if  D..  =  det  a^. 

and  if  all  that  was  known  about  the  elements  a—  was  that      ja.  .  |     ^     M  for  all  i,  j 

with  M  a  positive  constant,    then  the  most  one  could  state  is  the  result  of  Hada- 
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mard'  s  Theorem,       |DN  |  ^  N2     M     .      This  implies  that  |DN|  ^  MN2    — >-  oQ 


6.      Cf.    Whittaker  and  Watson,    Modern  Analysis ,    Cambridge  University  Press, 
Fourth  Edition,    1940,    p.    212. 


as  N  ->-   oO       •     Thus,    we  would  have  no  guarantee  that  a  thermodynamic  free 
energy  per  particle  would  exist  independent  of  N,    and  Klein'  s  Theorem  could  not 
be  applied. 

A  condition  on  the  elements  a- •  suf/Eicient  to  insure  the  existence  of  the 

required  limit  is  easily  found.     The  antecedent  of  Hadamard"  s  inequality  is  that 

7 
of  Kelvin,    Muir,    and  Hadamard     which  states  that 


7.     Hardy,    Littlewood,    and  Polya,    Inequalities,    Cambridge  University  Press, 
1934     p.    34 

Hence  if  2L    |a     |        ^    M     for  all  j,    then   |dn  |  ^  MN  so  that  DN  =  0(MN) 

fc=l        Jk 

if  M  exists.     In  our  case,  , 

M  ^  -2       2 

oo  J  2 

a   XZraDfeT)-2f4:(r)r2[i-e^[-Ah(r)i]cir    , 


where  f(r)  is  the  radial  distribution  function  of  ions  about  a  given  ion.     This  is 
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1/3  ? 

an  average  value  of  the  required  sum.     As  r-^  qQ    ,    h(r)  AJ  r/V  and  f(r)  /\»r    . 

Thus,    the  integral  does  not  converge  as  N,    V— >-0Q   ,    N/V  finite.     However,    we 

can  make  the  integral  converge  by  introducing  a  summability  factor  expj  -uh(r)  [ 

a  a  smallj  real,    positive  constant,    into  the  fjrce  law,    Eq.    (3).     Then 

Once  it  is  assumed  through  this  factor  that     £_|a  i I      is  bounded,    then  D^r  = 

I  IK'  " 

0(M    )  and  we  are  justified  in  employing  Eq.    (24)  for     \DN\V. 

Returning  to  Eq.    (19)  for  G     ,    we  now  have  the  result  that  the  singula- 
rity of  Gj^t  is  a  branch  point  of  order   1,    and  in  the  limit  N-^-oCi    occurs  at 

-1 

x  =   L(2DkT)      \      . 

v  '  m 


We  follow  Kramers'    method  for  determining  the  limiting  function 
1/3 

<*-'=- ixZDM"      ,    <DN>.=.(-ix)  AN    . 


Q  (x,    ni/J/T).     Let  us  put 


(28) 


Therefore, 


A,  =  Al^  +  ^'W    *   hi~-°  ■ 


(29) 


Note  that  A-.  depends  on  x  only  through    oC, 

Differentiating  AN  with  respect  to  o(  ,    and  noting  that  the  spatial  aver 
ages  of  the  minors  of  the  diagonal  elements  of  A^,  refer  to  a  system  of  N  -   1 
ions  and  that  these  averages  are  all  equal,    we  have 


< 


IA«    =  na      -  NAM_,     •  (30 
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Since    o(_^>  T      ,    and  since  Klein'  s  Theorem  shows  that  A„  depends  only  on  the 

1  /3 
variables  N  and  n    '     /T,    let  us  write  .1 

where 

y  =  n</  _  n  =  N/V,  X        =  an1/3, 

j  m 

with  a  a  numerical  constant.     Consequently,    we  are  to  find  the  function  0}(y) 
which  satisfies  the  differential-difference  equation  (30)  for  N  sufficiently  large. 
Using  the  form  for  A      in  Eq.    (30),    then 


where  y>  =   f,1^  =  [^-l)/\/]  ^     «     ^     ~     ^ /N)     ' 

However, 


Therefore, 

<34) 

/a  '/3  1/3  / 

l_  a(y'/3  -  I  -a.y      +    <llj/3/3N     ---     . 

Thus,    using  Eq.    (34),    we  have  in  the  limit  N  — >-  OO    ? 

3y^    =oJ    -.xfrfcd/.o)^]    ,         (35) 

-5 

when  |y  |     <^    1 /a    . 
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The  exact  solution  of  Eq.    (35)  can  be  conveniently  represented  in  para- 
metric form.      Defining 

F(y)=     -L^/OJ)     ^        > 
we  have 

Uj=e7('  +  3F)  >      k^'F'Ai  +  ^f)3,        06) 

where  K  is  an  integration  constant,     Kramers  has  neatly  determined  this  con- 
stant to  be  , 

K   =  i  C4ir)'  07) 

by  finding  A-^,  for  small  concentrations.  It  could,  of  course,  be  found  by  fitting 
the  results  to  the  Debye-Hueckel  limiting  law.  Eq.  (36)  differs  from  Kramers' 
result  in  that  y,    through  ^C,    is  a  function  of  x. 

Consequently,    through  Eq,1  s  (36),    (31),    (28)  and  (24),    we  find  that 

exf>(£)   =  (-tx)£F/(l+3F)    .  08) 

This  then  yields 

>/'  F  i^L 

G    =[|-(-iX21>bT)Xl      &  +  3r;e     /(-UU         ,  (3,) 


where  \        is  to  be  found  from  the  corresponding  singularity  of    Q 
Setting 

Z=-Lxzr^         }      <r%  K  n  (rj/PkT)       } 


(40) 
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F  as  a  function  of  z  is  determined  by 

fV(i+sf)*  -  *J/*3 . 


(41) 


(42) 


Using  the  value  of  A  in  Eq.    (10)  for  N  large,    Eq.    (17)  becomes 
i-l  00 


where  Eq ,    (39)  has  been  used  and  the  line  z     =  constant,    z      real  and  positive, 
is  to  the  right  of  the  singularities  of  the  integrand. 

The  integral  I  is  evaluated  by  the  method  of  steepest  descent.     The  de- 
tails are  presented  in  the  Appendix  and  we  shall  use  the  result  here. 

It  turns  out  that  I  has  two  forms  depending  on  the  value  of     (J    .     For  N 
sufficiently  large,    we  obtain 

-A\p/HkT  =  rzz  ^  o^<rc  ,    (43) 

-A^hT-|(r/0-iln^O-f      ,    o*c   , 

(44) 
2/3 

with     Yc  =  (2/3)      . 
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Discussion 


We  are  now  able,    through  Eq,    (43)  and  (44),    to  calculate  all  of  the  ther- 
modynamic properties  of  the  mixture  of  ions.     However,    we  will  restrict  ourselves 
to  the  more  important  quantities.     It  is  to  be  recalled  that    ^  =  (4Tf  f?)  Y*>  yl  I  U  R  '   J 
where    f\     is  related  to  the  ionic  strength  through  Eq.    (6). 

For  the  decrease      ^P  =  P  -  P  in  the  osmotic  pressure  P  we  have: 


(kT/2i/)y 


3/Z 


x<K 


Ap,-1L(^/H)=   I    -(kT/£^/0-t]     ^*Xt 


(45) 


where  v  =    1/n.      Since  P  =        (kT/v),    the  osmotic  coefficient  is  given  by 

r  3/Z 


1  -g    = 


,«*L 


i[(tfO-i]       ,y*rc 


(46) 


A  plot  of  the  osmotic  coefficient  against  the  square  root  of  the  concentration  is 
shown  in  Fig.    1.      We  obtain  the  Debye-Hueckel  limiting  law  up  to    the  equivalent 
concentration    Q        .      Deviations  from  the  limiting  law  in  the  proper  direction  set 
in  for  concentrations  greater  than  the  critical  concentratbn. 

For  the  decrease     £\  (D      in  the  thermodynamical  potential  we  have: 


A(p/N  =  ^/N+^p=i  kJ[zm.ium)-c\  ^  «« 


The  logarithm  of  the  activity  coefficient,    fa^,    for  an  ion  of  the  ith 
species  is  defined  by 
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l»fai    '0/feT)|^      '  (48) 


We  then  find: 


'»  "fki  - 


mine< 


The  nature  of  the    singularity      in  Z\U>*    at        Q    -  ^tf^  is  easily  deter- 

:d  to  be  a  discontinuity  in    (  O  A^f  /otf    /v-"      .     This  would  mean  that  the  osmotic 

C 

coefficient  is  continuous  with  a  continuous  first  derivative  at     q^     .     It  is  the  spe- 
cific heat,    however,    which  exhibits  the  most  unusual  behavior. 


We  easily  calculate 


(50) 


Defining  the  change  in  specific  heat  by 


AC    =   ^j-<>E)      , 


we  find: 


C3k/4)YSA       ,  V4*c  (S1, 


)     »  —  "C 


Consequently,    the  modified  theory  yields  a  constant  change  in  specific  heat 
above     fi        with  a  discontinuity  in  the  slope  at     0r    . 
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It  is  to  be  emphasized  that  although  the  expressions  for  the  common 
thermodynamic  functions  change  their  analytical  form  at        o     —   vr>  >    they 

are  still  continuous  through  the  critical  concentration. 

The  fact  that  the  modified  theory  yields  a  limiting  partition  function  for 
all  concentrations  of  the  ions  as  the  ionic  radius  approaches  zero,    demands  that 
a  revision,    and  in  some  aspects  a  reversal,    of  Kramers'    conclusions  be  made. 
The  primary  assumptions  underlying  the  theory  are:     1)  point  ions;  and  2)  small 
fluctuations.     Although  the  same  assumptions  underly  Kramers'    theory,    the  as- 
sumption of  small  fluctuations  is  not  adequately  treated  in  his  theory  because  of  the 
use  of  the  Gaussian  charge  distribution.      The  selection  in  the  modified  theory  of 
those  systems  in  Kramers'    ensemble  having  a  given  ionic  strength  is  consistent 
with  the  assumption  of  small  fluctuations. 

We  believe  that  the  existence  of  the  discontinuity  in  the  modified  theory 
is  due  to  the  "association"  of  ions,    although  we  are  unable  to  demonstrate  this  by 
a  pertinent  calculation  (except  through  the  behavior  of  the  osmotic  coefficient,    for 
example).     It  is  then  implied  that  fluctuations  become  significant  for      q     >  q\    , 
However,    the  association  is  not  the  catastrophe  of  the  Kramers'    theory.     We  may 
then  conclude  that  the  thermodynamic  functions  exist  in  the  limit  of  zero  ionic  ra- 
dius at  concentrations  greater  than  the  critical  concentration.     A  further  implica- 
tion is  that  when  fluctuations  become  important  it  is  also  required  to  bring  into 
the  theory  additional  properties  inherent  in  the  existence  of  discrete  charges,    i.e. 

that        y^    e**z.      =  constant 

Our  interpretation  of  the  modified  theory  is  that  for  concentrations  less 
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than  the  critical  concentration  fluctuations  are  unimportant  and  that  deviations  from 
the  limiting  law  are  indicative  of  the  effect  of  finite  ionic  size  and  of  the  approxima- 
tion of  the  ion-dipole  (solute-solvent)  interactions  by  a  macroscopic  dielectric  con- 
stant.    For  Concentrations  greater  than  the  critical,    association  plays  an  important 
part  and  this  involves  not  only  the  finite  ionic  size  and  dielectric  constant  approxi- 
mation but  also  the  fluctuations,      The  modified  theory  does  show  the  proper  direc- 
tion for  the  deviation  of  the  osmotic  coefficient  from  the  Debye-Hueckel  law  for 
concentrations  greater  than  the  critical, 

The  critical  concentration  of  the  modified  Kramers     theory  is  27  times 
greater  than  the  critical  concentration    of  the  Kramers'    theory  (Kramers      results 
correspond  to  zS  l),      For  a   1-1   salt,    the  new  critical  concentration  is  approxi- 
mately 0.81  mols/liter   rather  than  0.03  mols/liter.     The  modified    theory  does  not, 
therefore,    allow  us  to  conclude  that  the  association  of  point  ions  becomes  import- 
ant at  small  concentrations, 

The  result  of  the  work  of  Bjerrum  and  Fuoss  that  association  does  be- 
come important  at  small  concentrations  is,    of  course,    in  qualitative  agreement 
with  the  observations.      The  present  theory  makes  it  clear  that  one  cannot  avoid 
the  introduction  of  the  finite  ionic  size  if  even  the  qualitative  behavio  r  of  the  small 
concentration   deviations  from  the  Debye-Hueckel  law  is  to  be  reproduced       The 
discontinuity  at      Q         probably  does  not  actually  exist  in  that  a  finite  ionic  size 
plus  fluctuations  undoubtedly  force  association  to  be  gradual  rather  than  sudden, 

even  though  the  suddenness  is  not  catastrophic. 

8 
Recently,    J    E     Mayer      developed  a  statistical  mechanical  theory  of 


1.      J.E.    Mayer,    J.    Chem     Phys.    18,    1426(1950). 
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strong  electrolytes  which  is  based  on  the  Ursell-Mayer  cluster  integral  expansion. 
His  approach  is  probably  more  suitable  than  ours  if  one  wishes  to  take  the  detailed 
ion-dipole  interaction  and  van  der  Waals  forces  between  ions  into  account.     Since 
it  is  practically  impossible  to  evaluate  all  the  required  cluster  integrals,    Mayer 
limits  his  main  calculation  to  finding  the  contribution  of  the  so-called  ring  clusters 
(in  which  each  particle  interacts  with  only  two  others)  to  various  thermodynamic 
quantities.      This  contribution  seems  to  be  the  part  of  the  complete  expansion  which 
gives  rise  to  the  Debye-Hueckel  approximation.     Since  we  obtained  similar  results, 
it  appears  that  sphericalization  (the  condition  expressed  by  Eq ,    (7))  and  the  neg- 
lect of  all  cluster  integrals  except  rings  give  almost  equivalent  approximations  t_ 
the  partition  function. 

It  has  been  emphasized  by  Mayer  that  the  long  range  coulomb  forces  lead 
to  difficulties  in  the  evaluation  of  the  partition  function  of  a  system  of  ions.     The 
difficulties  are  avoided  in  his  theory  by  the  replacement  of  the  coulomb  law  with 
(£.   Gj   /r)  exp(-ur)  where  p.  is  a  very  small  positive  number,     At  short  distances 
this  expression  reduces  to  the  coulomb  law,    but  at  large  distances  where  ele 
cal  neutrality  prevents  a  physical  divergence  in  thermodynamic  quantities,    the 
factor  exp  (-|J.r)  prevents  a  mathematical  divergence  in  the  partition  function-. 
Mayer'  s  final  results  become  independent  of  u  as  u  — -^  0 ,     The  present  analysis  has 
also  shown  the  utility  of  the  convergence  or  summability  factor  exp  (-ur).     It  was 
useful  not  only  to  prevent  a  possible  divergence  but  to  insure  the  existence  of  a 
limiting  free  energy  per  particle  which  is  necessary  for  the  application  of  Klein'  s 
Theorem.     As  in  the  Mayer  theory,    if  one  takes  into  account  the  electrical  neu- 
trality at  every  step  in  the  analysis   it  should  not  be  necessary  to  introduce  tne 
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convergence  factor.     However,    since  this  requires  a  much  moi\e,  extensive  analysis 
we  have  avoided  it. 

The  above  analysis  had  its   origin  in  the  work  of  M,    Kac  on  the  Gaussian- 
model  of  a  ferromagnet  and  of  M.    Kac  and  the  authors  on  the  spherical  model  of 
cooperative  phenomena. 


9.     Kac        M.    and  Berlin,    T.H.,    Phys.    Rev,    75      1298(1949).     Montroll,    E 
Nuovo  Cimento  Supplemento  VI,    266  (1949). 


The  authors  wish  to  thank  Mr.    L.    Witten  for  several  discussions  and  his 
aid  in  preparing  the  Figures. 


APPENDIX 


We  shall  present  here  the  details  of  the  evaluation  by  the  method  of 
steepest  descent  of  the  integral  in  Eq,    (4>2^„ 

V'00  ri/2r  nN/2 

J  a(2irt)-'  jd^[/-(Am^/DkT^]        [r'(i+3F)exK-F+l)J 


?o-LC0 


(A   1) 


where  F2/(l  +  3F)3  =     Y  /Z3, 


and  the  line  z     =  constant,    z     real  and  positive,    is  to  the   right  of  the  singulari- 
ties of  the  integrand. 

The  function  F  is  a  three-valued  function  of  z.     We  can  construct  a  Rie 

mann  surface  of  three  sheets  to  make  F  an  analytic  function  of  z. 

3 
Let   ^-(2/3%)  •       The  solutions  of  the  cubic  equation  for  F, 


are: 

2,l/3  ..  l/3 

F    (z)=  ^(f-0    ~U)     L        -   U3M  ' 

F2<z>  =  s(f-0  -60  L.'3-oo2M  3     5  (A  2) 

the  subscripts  denoting  the  different  sheets  of  the  Riemann  surface,    and  where 

oj  =  e.*\.{LMlt)    , 

L  .-(f/27)(:f-a?  +  |)+o(?/27X?-|)'^  (A3) 

The  discriminant  of  the  cubic  is       £  4  *?  /2  7/w         4  y    ' 
The  branch  points  occur  when  the  discriminant  equals  zero  and  are    ^  =  0  .  -j-  > 


A-2 


and  this  means  that  the  cubic  has  multiple  roots,      Note  that  the  point  at  infinity- 
is  a  branch  point  since  L— ^-M  as    \f\  -*■  00       ■     The  values  of  z  at  the  branch 

2.1  Z  2/3  */*      2 

points  are  z  =  0,      3X(3/2)  3  if  (3/2)      6J    _,  and     3^(3/2)      tU     • 

/     /    \2//j 

Consequently,    in  the  expression  (A  1)  for  I,    we  must  take  z     >  3^(3/2}        , 

Furthermore,    X        corresponds  to  this  largest,    real,    and  positive  singularity, 
and  we  conclude  that 


or 


m 


1/3     1/* 

3  -a     n     . 


By  constructing  a  cut  in  the  z-plane  as  shown  in  Fig,    A   1,    F .  is  ana 


J 


lytic  in  the  cut  plane, 


— 00-*- 

Z; 

1 

0 

*. 

Z4.=3y(3/2) 

3 

V 

Fig.   A  1.      The  branch  cuts  in  the  complex  z-plane 


Let  us  write 


G         -1.  .     -F+z 

e     =  z      (1   +   3F)e 


jr 


(A  5) 


G  =  z-lnz-F+  ln(3F  +   1), 
and  attempt  to  find  a  saddle  point.      We  have 


A-3 


dG/dZ        =  z_1(z  -   1    -  3F)  (A  6) 

and 

d2G/d*2      =  2(1  +  3F)2/z2(2  -  3F).  (A    7) 


Letting  z     be  such  that  (dG/dz)        =  0,    then 
s  s 

zs  =   1  +   3F(zs).  (A  8) 

Using  the  relation  between  F  and  z,    we  obtain 

3/2  2/2 

FK) =  tX     y    zs  =  l ±3*  <A  9> 

Also 

(d2G/dz2)        =  2/(3   -  z   ).  (A  10) 

zs  5 


We  conclude  that  a  saddle  point  exists  for    jf  in  the  range  0  ^  o   ^  00 

ddle 

/     2/3 

r     ^     in  the  range  0  <    )f  <    ^c  =  W3)  ■     When     )f  ss   % '     I,    z+  =  B, 


when  z     =z     =   1    -  3  o  >    and  that  a  saddle  point  exists  for  z     =  z.    =   1  +  3  0 

s  -  .      r  s 

fo 


(A  11) 


and  the  second  derivative  of  G  becomes  infinite.     For     0    —     0,    zs —  1,    and  since 

^  =  (H/jy)   j       lf|^(.3¥)"       ~  00  .  It  is  easily  found  that  when    |^|-»  00, 

F3C*J=  1  -1    "  ^?)"'-  "-  • 

Therefore,    the  branch  Fj(z)  corresponds  to  F(z    )  =  +  Q  and  the 

branch  F^(z)  corresponds  to  F(z    )  =   -  q  . 

From  Eq.    (A  5) , 

G(zJ  =  zc   -  F(zJ   -   1  +  2F(zs), 
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so  that 

G(z.)  =   1  +2  ^         and  G(z_)  =    1    -  2  fl         •  (A   12) 


Since  for  large  N,    I  is  essentially  exp  I  NG(z    )/2  J    or  exp  InG(z_)/2J    ,    it  is 
plain  from  Eq.    (42)  of  the  text,    where 

y*    -H/2 

ex|>[-^/feT]  =  Otn)    e       I  , 

that  the  significant  contribution  to  I  arises  from  the  saddle  point  z ,  .     As  a  con- 
sequence    the  integral  I  is  to  be  evaluated  on  sheet  1  where  we  require  F(z)  S 
F    (z).     This  means  that  we  can  now  write 

*+  +  ^G  i  -i/2  2 

I  a  (ziuf  )<L*  [l  -  3  r(S/2)   z" ']  ex  j,  |  [S  CO*  ±  0&/d  i'^-  *)  -  -J 

so  that 

X  a [j --5/(3/2)%-']    frN(dtyd*%+]"e*f>[NGC2+)/2]  •  (A  ,„ 


Hence  f 


and  for  N  very  large, 

*A  2/3 

-  A^/MfeT  *  Y  >  ^  ^  -  (2/3)      .  (a  is) 

As        J    — ^     Of  >    the  saddle  point  z,-^  3,    a  branch  point  of  F,    and 

the  derivative  Id   G/dz^ J   -**  QQ     .     It  would  appear  that  the  branch  point  limits 
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the  validity  of  the  integral  I  because  if  we  use  the  above  saddle  point  equations 
for  Q     >     o^  ,    then  z    >    3,    the  second  derivative  of  G  changes  sign  and 

the  path  of  steepest  descent  is  along  the  real  axis,       This  would  make  I,    and  /\  XL) 
purely  imaginary  as  can  be  seen  by  analytically  continuing  Eq.    (A  14).     The 
apparent  limitation  is  equivalent  to  the  breakdown  found  by  Kramers.     His  for- 
mulas become  meaningless  when       Y  =   jl       and  the  branch  point  is  reached, 
because  his  equations  become  undefined  for     o    >  cTr    • 

The  path  along  the  real  axis  through  the  saddle  point  z     is  not  mathema- 
tically legitimate  for      q  >X,       because  the  imaginary  part  of  G(z)  for  this  path 
is  identically  zero.     This  path  of  steepest  descent  can  not  be  deformed  into  the 
original  path,    the  line  for  which  the  real  part  of  z  equals  z    .      We  therefore  sup- 
pose that      q    >   q        and  investigate  in  further  detail  the  behavior  of  the  integrand 
in  the  neighborhood  of  the  branch  point  of  F,   P     =  9/4. 

Let  us  set    v  =  (9/4)  +  €.   with    (  €.  \    very  small.      Then,    since  z  =  3  X  S 


3(r/^[i  +(4/*7)e ]  . 


Furthermore, 

.  i  /o 

P(i)  =(2/3)  -C4'3,/2/9)e'%  (40/8 f)€  -  (59Z^Z/z43)e    . 

From  the  definition  of  G(z),    we  then  obtain 


'       (A  1  6) 


where 


GO)  =  G0  +  G^g   +*e      -    -  >  (a  i?) 

G0=  3(*/rc)-  ln(^r/rc)   -  2/3  , 

Go-  C4/9)[«/arc)-i]  ,  Ul8) 
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Now  when     q    •>  Q       ,    we  note  that  G'      >    0.     Consequently,    for     |£| 
small  and  real  part  of   £,    negative,    G(z)  has  a  maximum  value,    G    ,    at 

€.   =   0,    z  =  Zj   =  3  (jf /  0^  )         •     The  remaining  question  concerns  the  possi  - 
bility  of  locating  a  path  of  steepest  descent,    imaginary  part  of  G(z)  =   0,    which 


i  6 
can  be  deformed  into  the  original  path.     If  we  set    £  =  re       with  r  small, 

3/2 
imaginary  part  of  G(z)  =  G'     r  Sin  6  +   o(_r    '     Sin(36/2). 


(A  19) 


Setting  this  equal  to  zero,    we  have  for  the  possible  path 


Sm6 


(A  20) 


The  path  cannot  cross  the  cut  from  0  to  z     =  3()f/^c_  )         in  Fig.    A   1.     There- 
fore,   in  the  upper  half  plane  we  take     Q  =."{[— Q>     and  in  the  lower  half  plane  we 
take      ©  =  —71"  +  CD  where  <3>    is  positive  and  small  in  both  cases.     Hence 


^--^»^§^)^G^}^ 


(A  21) 


and  a  path  not  on  the  real  axis  exists  on  both  sides  of  the  cut.     The  path  is  shown 


in  Fi  g.   A  2. 


-  00 


e=ir 


Q  =  -f 


-l^lane. 


2     C 


1®. 


^i 


Fig.   A  2.     Path  of  steepest  descent  for  J  ><L 
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In  the  neighborhood  of  Zj,    the  integrand  becomes  infinite  like 

Q-(3r/2rc);T'/2*  (zrAe/2 , 

and,    as  a  consequence,    the  contribution  to  I  from  an  infinitesimally  small  path 
circling  z,    is  zero  in  the  limit.     The  integral  may  then  be  written 


I 


,-J 


C2TT6f'^4)r/9OdH27/40^x(>[rM/OCG0H-GjG-Oj    .    (A22) 


c+c 


I  '  "-z 


For  N  very  large  the  maximum  is  very  sharp  and  the  contributions  to  I 
arise  from  the  parts  of  C.    and  C?  very  near  the  branch  point,     Therefore,    the 
contributing  parts  of  Cj  and  C^  may  be  replaced  by  integrals  along  the  real 

-tlT 

axis .     In  C,   we  may  use     ^  =  )r£  qq   -^    f  >.  Q    ,    and  in  C,  we  may  use 

e~reLir^  o^r  <  °°-  Hence0 


I  a  mfX^KW^^'^     exH(N/2)(Se-^r)j 


(A  23) 


+  (ztLT1 


OQ 

(4y/9rc)(27/4)'/2{(-*)r"'/Vlff Ze^]Hli)(G0  -Q/)]  , 


and 


,^l* 


I 


^^(2/3/V/O^O"'  e^[HGj2]  . 


(A  24) 
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:|=[-^AT] *  ^(3/z)/V/irt)L0f/)rc)-iJe)()=[(N/2)(Go- ')]  >    (A  25) 


and  for  N  very  large, 


We  note  from  Eq.    (A  15)  and  (A  26)  that    /4  sP      is  continuous  when 


-ApMxr  =2/3       ,    y=^c  • 


(A  27) 


